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Abstract 

We discuss the influence of the cosmological constant on the gravitomagnetic clock effect 
and the gravitational time delay of electromagnetic rays. Moreover, we consider the rela- 
tive motion of a binary system to linear order in the cosmological constant A. The general 
expression for the effect of A on pericenter precession is given for arbitrary orbital eccentricity. 
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1. Introduction 

Current cosmological models that take the acceleration of the universe into account involve 
a cosmological constant A of magnitude A « Hq/c 2 , where H is the Hubble parameter. 
Taking H ~ 70 km sec -1 Mpc -1 , we find that A ps 10 _56 cm~ 2 . The local effects of such a 
cosmological constant are expected to be very small. Nevertheless, great strides are being 
made in observational techniques in astronomy and it appears worthwhile to investigate some 
of the observational consequences of the existence of a cosmological constant. 

The dependence of the gravitomagnetic clock effect on the cosmological constant is exam- 
ined in section 2. Section 3 is devoted to the gravitational time delay and its dependence on 
A. The influence of the cosmological constant on the motion of a binary system is investi- 
gated in detail in sections 4 and 5. A bound orbit is shown to precess with a frequency given 
by |(Ac 2 / 'u)y/l — e 2 , where to and e are the orbital frequency and eccentricity, respectively. 
Section 6 contains a brief discussion of our results. 



2. Clock effects in the Kerr-de Sitter spacetime 



The gravitomagnetic clock effect is caused by the net spin of a gravitational source [1-4] . In 
general, the rotation of a massive body introduces a difference in the co- and counter-rotating 
orbital periods of test masses in motion around the central source. We wish to investigate 
this gravitomagnetic effect in the equatorial plane of the Kerr-de Sitter spacetime to find 
what influence the cosmological constant would have on the clock effect. The Kerr-de Sitter 
metric in the standard Boyer-Lindquist coordinates = (t,r,9,<f>) is given by [5, 6] 
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where units are chosen such that G — c — 1, unless specified otherwise, and 

E = r 2 + a 2 cos 2 9 , y = 1 + —a 2 cos 2 6> , 
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A = r 2 - 2Mr + a 2 - ^r 2 (r 2 + a 2 ) 



(1) 

(2) 
(3) 



The mass of the source is M, J = Ma is its angular momentum and A is the cosmological 
constant. The metric (1) represents a rotating Kottler spacetime [7]. The Kottler spacetime 
is also known as the Schwarzschild-de Sitter spacetime. Useful background information on 
de Sitter spacetime and the cosmological constant can be found in [8] and [9,10], respectively. 

It is interesting to note that in the Newtonian limit the physical content of equation (1) 
reduces to a Newtonian gravitational potential & N = —GMr' 1 — |Ac 2 r 2 , which satisfies 
the Poisson equation V 2 <3>at = 4irGp N with p N = M8(r) + p A . Here pa = — c 2 A/(47rG) is 
the effective negative uniform density of "matter" represented by the cosmological constant. 
The effective repulsive force due to this source is then given by |Ac 2 r per unit mass. On the 
other hand, the relativistic interpretation of the cosmological constant in terms of vacuum 
energy involves a "substance" with equation of state pc 2 + p = such that the density p is 
positive, but the pressure is negative and is given by p = —c 4 A/(8ttG). 

Let us now consider circular geodesies of this spacetime. The stability of circular geodesic 
orbits in the equatorial plane of Kerr-de Sitter geometry has been investigated by Howes 
[11]. It turns out that stable circular orbits exist out to a radius of (3M/A) 1 / 3 . On the other 
hand, the inner stability limits vary depending on the sense of orbital motion as well as the 
values of a and A [11]. For A = 0, we recall that for the corotating case, the inner stability 
threshold varies from r = 6M to r = M as a : — > M, whereas it varies from r = 6M to 
r = 9M in the counterrotating case. 

The geodesic equation for the radial coordinate reduces to [12, 13] 

^(|) + 2 ^(|)+^ = ° > ( 4 ) 

which has the solution 

- = a±u, K . (5) 
Here ojk is the modified Kepler frequency 

^ = -T-3 A > (6) 

which can be interpreted as being due to a net mass of M+ (47r/3)/)Ar 3 , where p\ = —A/ An 
is the constant Newtonian density mentioned before. In connection with equation (6) as 



well as other modifications of Kepler's third law to include dark matter, etc., we note that 
the motion of planets in the solar system can be used to study deviations from Keplerian 
motion; however, for A lCT 56 cm~ 2 the corresponding effects are too small to be detectable 
at present. 

There are three distinct clock effects that can be studied using counterrevolving circular 
geodesic orbits. To find the influence of the cosmological constant on the clock effects, 
we consider these in turn. Integration of equation (5) over 2n for corotating and — 2n for 
counterrotating test particles implies that t± = 2V±27ra, where Tk = 2-k/uok is the Keplerian 
period. Thus t + — t_ = 4na, just as in the Kerr case. Consider now an (accelerated) standard 
clock at rest on the circle of radius r in the equatorial plane. A standard clock measures the 
proper time along its world line. According to this clock, the difference in the time r' that 
it would take two free counterrevolving test particles to complete the orbit is 



<-7i = 47rayi-^-|(r« + a2) , (7) 

which is the observer-dependent single-clock clock effect according to the terminology of 
[12]. For the cosmological constant to have any significant influence on this clock effect, the 
orbital radius would have to be unreasonably large. In fact, A as 10~ 56 cm -2 has a totally 
negligible effect even at a planetary radius of r « 10 AU, since Ar 2 /3 ~ 10~ 28 . 

Consider now the proper periods of revolution r± of the two free standard clocks coun- 
terrevolving on the circular orbit of radius r; we find from equations (1) and (5) that 



r± = T K Jl- ^ - ^a 2 ± 2au K . (8) 



Expanding this expression in powers of e = a/M, we find that 

Ma 2 (l - Ar 2 ^ 



4na 

r + - r_ 



(9) 



2r(r - 3M) 2 

Thus the cosmological constant does not affect the observer-dependent two-clock clock effect 
to first order in a/M. 

Imagine now two standard clocks starting at the event characterized by (0, r, tt/2, 0), i.e. 

t — at (f) — 0, and moving freely in opposite directions on the circle of radius r. Since the 

corotating clock is slower, their first meeting occurs at <f>i = n — a, where a = clujk- Let 
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their nth meeting point be the event characterized by t n and <p n , where t n = jnT K (l — a 2 ) 
and 4> n = mr(l — a) modulo 2ir. The observer-independent two-clock clock effect refers to 
the difference in the proper times r+ — t~ of the clocks at their nth meeting, where 



, 1 . . / 3M A „ 
T ± = -nT K {l T a)^l- — --a*±2a . (10) 

To first order in a/M, this clock effect is also independent of A and is given by 

^r(r - 3M) 

The diametrical line joining the points of encounter with the origin of the circular orbit 
undergoes a precession in the opposite sense as the rotation of the source. The preces- 
sion frequency is given approximately by nna/r^, which to lowest order coincides with the 
precession frequency of a fixed torque-free test gyroscope in the equatorial plane at radius r. 

The gravitomagnetic clock effect (9) far from the source, 3M << r << (3M/A) 1 / 3 , is 
essentially independent of A, and is basically proportional only to a = J/M, which is the 
specific angular momentum of the source. That is, the effect has a topological character as 
it is practically independent of the radius of the orbit. It is also essentially independent of 
the gravitational constant G; therefore, the effect can be "large". For instance, for space- 
borne clocks around the Earth the effect is ~ 10~ 7 sec; however, the detection of the effect 
requires knowledge of satellite orbits to millimeter accuracy [3] . This seems to be somewhat 
beyond present capabilities by about an order of magnitude; moreover, the influence of the 
cosmological constant on the clock effect is too small to be detectable in the foreseeable 
future. Another aspect of the clock effect is that it takes longer to complete a prograde 
circular orbit in the equatorial plane than the corresponding retrograde orbit. Thus free 
motion in the same sense as the rotation of the source is slower than motion in the opposite 
sense, which is contrary to what one might expect from the "dragging of space" by the 
rotating source. 

There is a close relationship between the clock effect and circular holonomy. In fact, 
circular holonomy in the Kerr-de Sitter spacetime can be developed in close analogy with 
that in the Kerr spacetime [14]. This would involve the parallel transport of a vector around 
a constant-time circular 0-loop in the equatorial plane of the Kerr-de Sitter spacetime. Let 



us just mention here the result that there is a band of holonomy invariance for r* < r < 
(3M/A) 1//3 given by nf(r) = m, where n and rh are positive integers and 

2M 2Ma 2 M 2 a 2 1 , 1, ,/ 2M\ 

A ^ i A ^ f 1 i 1 



f\r) = 1 — - -Ar 2 + -Aa 2 1 + — (12) 

r r d r 4 3 3 V r / 

to first order in the cosmo logical constant. Here r* is the unique positive root of /(r) = 
for sufficiently small A. 

3. Time delay 

Let us consider the time delay in the propagation of electromagnetic rays in the Kerr-de 
Sitter spacetime. The influence of the cosmological constant on the local propagation effects 
is expected to be very small; therefore, it is sufficient to consider the linearized Kerr-de Sitter 
metric in isotropic coordinates on the background Minkowski spacetime. In this case, the 
time that it takes for the rays to travel from the point Pi : (t, pi) to the point Pi : (t, p?) 
can be written as 

*2 - h = \p 2 - Pi\ + Age + A GM + A A , (13) 

where p is the position vector in the background Minkowski spacetime and A GE and A GM 
respectively represent the standard gravitoelectric and gravitomagnetic time delays [15]. 
Moreover, Aa is the additional time delay due to the cosmological constant. This delay can 
be calculated using the approach developed in [15]; that is, 

A = Yc Q h °e kakf>dl > ( 14 ) 

where h a p = g a/ 3 — r\ a p , k a — (1, k), k is a unit vector along the direction of propagation 
of the ray in the background Minkowski spactime and dl is the element of straight line 
connecting P 1 to P 2 - 

Concentrating on the cosmological constant, we note that for M = a = 0, equation (1) 
reduces to the de Sitter metric. Under the coordinate transformation 

ri » - TTW ' <15) 

the de Sitter metric takes the isotropic form 
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where tp = Ap 2 /12. It follows from equation (16) that for the cosmological constant alone, 

^oo = ^Ap 2 , hm = , hij = -^Ap 2 5ij . (17) 
Using these results in equation (14), we find 

4 *4/>' (i8) 

It is straightforward to compute this integral and the result is 

Aa = ^T l^ 1 ~ P 2 \ (Pi + Pi ' P2 + Pi) • ( 19 ) 
36c 

We note that for A 10~ 56 cm~ 2 , A A ~ 2 sec for a ray of radiation crossing the disk of our 
Galaxy. 

The calculation of first-order gravitational delays in equation (13) assumes that the ray 
follows a straight line. Therefore, there must be in addition a geometrical delay that takes the 
actual deflected path into account; for weak lensing, the geometrical delay is of second order 
in the deflecting potentials. It is important to recognize that the cosmological constant does 
not participate in the deflection of light rays in the Kottler spacetime [16, 17]; therefore, the 
cosmological constant does not contribute to the geometrical delay in the Schwarzschild-de 
Sitter spacetime. 

It follows from equation (19) that locally, i.e. for distances much smaller than (3M/A) 1 / 3 , 
the time delay due to the cosmological constant is very small compared to the corresponding 
Shapiro time delay. On the other hand, equation (19) indicates that for distances approaching 
A -1 / 2 , the time delay could be rather significant, but then our linear approximation scheme 
may break down. Thus for cosmological observations Aa should in general be taken into 
account, since the present uncertainty in the measurement of gravitational lensing time delay 
is about | day [15]. 

4. Orbital precession 

The solution of the equations of motion of a test particle in the Kerr-de Sitter spacetime 
is quite complicated. On the other hand, it turns out that one can obtain the main results 



regarding the gravitomagnetic clock effect, etc., from the linear approximation of general 
relativity. Therefore, we consider a linear post-Newtonian approach using a Lagrangian 
of the form C = —m ds/dt, where m is the mass of the test particle and the metric is 
expressed in (t, p) coordinates, where p = (x, y, z) indicates a point in space with isotropic 
(post-Newtonian) Cartesian coordinates. The equations of motion then take the form 

d 2 P GMp 

-d¥ + — - F > (20) 

where F is the post-Newtonian perturbing function. The details of this process of reduction of 
the equations of motion to the form (20) are straightforward and have been described in detail 
in [4, 18]. Under certain circumstances, the orbital perturbations can be easily characterized 
in terms of a modified orbit as has been done for the gravitomagnetic clock effect [4]. To 
illustrate this procedure, it is instructive to digress here and consider the case of bound 
motion in the Kerr-Taub-NUT spacetime; the results can be compared and contrasted with 
those of the Kerr-de Sitter spacetime. Linearizing this metric in the angular momentum and 
the NUT parameter of the source, we find the Schwarzschild metric together with the Lense- 
Thirring and Taub-NUT off-diagonal terms. Introducing the isotropic Schwarzschild radial 
coordinate p and the corresponding Cartesian coordinates, we find that this gravitoelectric 
background is perturbed by a gravitomagnetic field, 



B, 



GJ_ 

cp 5 



3(p-3)p-p 2 3]-^f (21) 



due to the angular momentum (i.e. the gravitomagnetic dipole moment J) and the grav- 
itomagnetic monopole moment (—c 2 £) of the source. Here the NUT parameter t has the 
dimension of length. The post-Newtonian gravitomagnetic orbital perturbation can then be 
obtained from equation (20), where F is of the form F = — 2v x B g /c. The linear perturba- 
tion of the orbit due to the angular momentum of the source has been worked out in detail 
in [4] in connection with the gravitomagnetic clock effect. Therefore, we only indicate here 
the orbital perturbations due to F^ = —2c£p x v/p 3 using the method developed in [4]. 

For the motion of a test mass around a gravitomagnetic monopole, the energy of the 
particle is conserved and so is its angular momentum, which is however augmented by a 
contribution from the interaction of the particle with the monopole, i.e. j = mpxv—2mc£p/ p 



is the conserved quantity. Following the analysis given in [4] , let us note that the unperturbed 
orbit is given by a Keplerian ellipse in the (X, F)-plane of the (X, Y, Z) coordinate system 
that is related to the (x, y, z) system by a rotation, 

x = cos fiX — sinficosi Y + sinfisini Z , (22) 
y = sinfi X + cos fi cos 2 Y — cosfisini Z , (23) 
z = sin? Y + cos? Z , (24) 

where i is the inclination angle and Q is the longitude of the line of the ascending nodes. 
The Keplerian ellipse is then given by X = pcostp, Y = psiny? and Z — 0, where 

a(l — e 2 ) dip L 



1 + ecos(y? — g) dt p 



2 



(25) 



and L = ^GMa(l - e 2 ) is the specific orbital angular momentum of the unperturbed 
orbit. Here a, e and g are respectively the semimajor axis, eccentricity and argument of the 
pericenter of the unperturbed elliptical orbit. 

Once the gravitomagnetic monopole moment of the source is turned on at t = 0, when 
observations begin, the only change that occurs in the orbit to linear order is that Z is no 
longer zero but is given by 

z = \ u i , (26) 

c 1 + ecos(<p — g) 

where \x = —c 2 £/ (GM) is the dimensionless strength of the monopole and ip — ip at t — 0. A 
detailed examination shows that there is no gravitomagnetic clock effect in this case to first 
order in £. In this linear order, the net angular momentum j turns out to have componets 
jx = —2mc£ cos ip , jy = —2mc£sm(p and j z = mL . 

The gravitomagnetic monopole moment of the source presumably exists for all time; 
therefore, the procedure described above has to be re-interpreted in terms of an osculating 
ellipse. Namely, the position and velocity of the test mass at any given instant of time define 
an elliptical orbit that is momentarily tangent to the actual perturbed orbit. Therefore, the 
perturbed motion may be described in terms of the evolution of the osculating ellipse as in 
the Lagrange planetary equations. In the case under consideration here the observations are 
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assumed to begin at t — 0; hence, the elliptical orbit (25) is simply the osculating ellipse of 
the perturbed orbit at t — 0. 

Let us now return to the Kerr-de Sitter metric and linearize it in the small parameters 
e = a/M and 5 = AM 2 with M^O. We then need to introduce an isotropic radial coordinate 
p that puts the Kottler metric in isotropic form. The connection between r and p is given 
by the differential equation 



dr\ 2 r 2 / 2M A 2 

1 r 2 . (27) 



\dp ) p 2 
To first order in 8, we find that 



2 



*>-' K -H-S '(2) ' 



where I(x) is given by 

«•>-/" 5 ■ w 

This integral can be evaluated and the result is 

128 16 1 

2J(ar) = K + x 2 + 16x + 601nx -, (30) 

x — 1 X X 2 

where if is a dimensionless integration constant that we can set equal to zero for the sake of 
simplicity. A detailed examination of the resulting equation of motion then shows that the 
dominant effect of the cosmological constant would simply be due to the effective Newtonian 
perturbation of the form F = Ap, where A = Ac 2 /3. Under such a perturbation the orbital 
angular momentum is conserved and the orbit is thus planar. On the other hand, the 
orbital energy is augmented by a contribution from the cosmological constant so that the 
conserved quantity is E, where 2E = v 2 — Xp 2 — 2GM/ p. Though the perturbing function is 
simply proportional to p in this case, it turns out that the simple method developed in [4] 
is not directly applicable in this case for orbits of arbitrary eccentricity, since the perturbing 
function has the form of an infinite series in powers of the eccentricity (cf. section 5). This 
is in contrast with the Kerr-Taub-NUT case. It is therefore necessary to adopt a different 
approach to the motion of the perturbed orbit unless e << 1, in which case the method of 
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[4] is adequate. The general solution to first order in the cosmological constant is considered 
in the next section; in the following, we study the average motion of the perturbed orbit. 

It is interesting to compute the average rate of precession of the pericenter under the 
influence of A. To this end, let us introduce the Runge-Lenz vector of the perturbed orbit 

Q = v x L - GMp/p , (31) 

where L = L Z and we note that for F = 0, Q = GMe(cosg, sing, 0) in the (X,Y,Z) 
coordinate system. For F = Xp, the orbit remains in the (X, F)-plane and Qx = XL Y, Qy = 
—\L X and Qz = 0. Thus in this case 

^ = -AL x p . (32) 

Averaging this relation over the "fast" motion, i.e. the unperturbed orbital motion of the 
test particle (with period T) would reveal the "slow" (secular) behavior of the Runge-Lenz 
vector over a long period of time. Defining the average over a period as 

T Jo 2it Jo 1 + e cos (ip — q r 



we obtain 



3 3 
< X >= — -aecosg , < Y >= — -aesing , (34) 



where we have used the relation 

1 f(o+2n cos 



271-4, (1 + ecosC) 3 2 (1 -e 2 ) 5 / 2 

Thus we find that 



(35) 



<^>=nxQ 0) (36) 

where 

2 GM 2 uj K ' 

Here u = 2tc/T is the Keplerian frequency of the unperturbed orbit. For A lCT 56 cm -2 , 

the precession frequency (37) is too small to be observable at present in the solar system 

[16]. The precession angle after one period is given by 
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To first order in eccentricity, this result agrees with previous work [16, 19], while the exact 



dependence on eccentricity as yl — e 2 differs from (1 — e 2 ) 3 given in [19]. 

To gain a deeper understanding of the average behavior of the orbit, we consider the 
osculating ellipse in this case and note that the orbital elements of this ellipse vary with 
time according to the Lagrange planetary equations [20] as 



da 2Aaev / l — e 2 sm(ip — g) 

dt oo 1 + ecos(y? — g) ' 

de 1 — e 2 da 

~dt ~ 2ea ~dt ' ^ ' 

dg _ Ay / (1 - e 2 ) 3 CO s(y9 - g) 



dt ecu 1 + e cos(</? — g) 

Employing the averaging method once again, we find that 



(41) 



da de . . 

< * > =°- < * > =°- (42 > 
< <k >= | Ayr^ , ,43) 

at 2 

where relation (35) has been used. The rate of advance of the pericenter given by equation 
(43) is consistent with our previous result. Thus in the presence of a cosmological constant 
the orbital plane does not change and the orbit keeps its shape on the average while precessing 
with frequency Ac 2 Vl — e 2 /(2oo). 

It follows from equation (37) that the frequency of the pericenter precession caused by the 
cosmological constant is directly proportional to the period of the binary orbit, which must 
therefore be very large for the effect to be perceptible. The ratio of this precession frequency 
to the Einstein pericenter precession frequency is given by Aa 4 c 4 (l — e 2 ) 3 / 2 / (6G 2 M 2 ), which 
is ~ 1CT 14 for the motion of the Earth around the Sun. We conclude that the influence of the 
cosmological constant on orbital motion is too small to be measurable in the solar system. 

5. Perturbed Keplerian motion 
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It is interesting to give the general solution of equation (20) for F = Xp to first order in 
the cosmo logical constant. For a positive cosmological constant, the effective force on the 
test mass is along the radial direction away from the source. The orbit is thus planar and 
p 2 dip/dt = L is the constant specific orbital angular momentum. In terms of u — p^ 1 , the 
radial part of equation (20) reduces in this case to 

d 2 u _i A , AA . 

+ u = p - 44 

dtp 1 LqU 4 

where p = L\j (GM) is a positive constant. The right side of equation (44) can be expressed 
as a function of <p if we substitute for u the unperturbed solution in accordance with our 
linear perturbation scheme. The unperturbed solution is given by 

u = - (1 + ecosy?) , (45) 
p 

where we have chosen the planar Cartesian coordinates such that g = for the sake of 
simplicity. Here e > is the orbital eccentricity of the conic section. For an ellipse e < 1 
and p = a(l — e 2 ), where the pericenter distance from the focus is p m i n = o(l — e) and a is 
the semimajor axis. For a parabola e = 1 and p m i n = p/2 is the pericenter distance, while 
for a hyperbola e > 1 and p = a(e 2 — 1), where the pericenter distance from the focus is 
Pmin = o(e — 1). Substituting equation (45) for u in the right side of equation (44), the 
resulting equation can be transformed to the form 

~. d 2 U _dU . 

< 1 -Vw- ( « +v = wkr- (46) 

where £ = cosy?, U = u — p~ x and q = — Ap 3 /Lg. Equation (46) can be integrated once and 
we find 

where C is an integration constant. Equation (47) can be written in the form 

d^ { ^ J 2e(i-a 3/2 (i+e0 2 (i-a 3/2 1 ] 

and integrated using formulas 2.264, 2.266 and 2.268 given in [21]. The result for an elliptical 
orbit is 
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2(l-e 2 ) 



- 3(l + eXsin<^) 



(2e H — )cosip + C cos tp + S sin y9 , (49) 



1 + e cos y> 



e 



where S is a constant of integration and 



J = 



1 



arc sin 



( 



1 + e cos (/? 



e + cos cp 



) 



(50) 



The integration of equation (48) can be similarly carried out for parabolic and hyperbolic 
orbits using the relevant formulas of [20]. We expect that u reduces to the unperturbed 
solution uq as A — > 0; therefore, it is useful to set C = qC + e/p and S = qS'. The constants 
C and S' can be uniquely determined using the initial conditions, i.e. the initial position 
and velocity of the test mass. In this way, we find the orbital radius p = u^ 1 = {p^ 1 + U)^ 1 
and then ip(t) by integrating dip/dt = L p~ 2 . 

6. Discussion 

In this paper we have considered some of the local physical consequences of the existence 
of a cosmological constant. In particular, we have considered the motion of test bodies and 
the propagation of rays in the Kerr-de Sitter spacetime. While the clock effect is rather 
weakly affected by the presence of a cosmological constant, there is a (gravitoelectric) time 
delay in the propagation of rays directly proportional to A. Though this effect is too small 
to be measurable in the solar system, it should in general be taken into account in the 
interpretation of gravitational lensing time delay. Moreover, the average behavior of bounded 
orbits has been examined and the general rate of orbital precession has been determined for 
an arbitrary eccentric orbit. On the average, the orbital plane remains invariant and the 
orbit keeps its shape but precesses due to the presence of a cosmological constant; however, 
the precession rate is too small to be detectable in the solar system. 
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